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Abstract 

It is assumed that, for weak spacetime curvature, the main grav- 
itational effect of the quantum vacuum stress-energy corresponds to 
adding two terms to the Einstein-Hilbert action, proportional to the 
square of the curvature scalar and to the contraction of two Ricci 
tensors, respectively. It is shown that compatibility with terrestrial 
and solar systems observations implies that the square roots of the 
coefficients of these terms should be either a few millimeters or a few 
hundred meters. It is shown that the vacuum contribution increase 
the stability of massive white dwarfs. 

PACS numbers: 04.50.+h ; 04.25.Nx 

1 Introduction 

In the study of quantum fields in curved spacetime it has been stablished 
that the quantum vacuum gives rise to a finite, non-zero, energy [1], [2]. Fur- 
thermore, some effects have been attributed to the gravity of the quantum 
vacuum, like the observed acceleration in the expansion of the universe [3]. 

The gravitational effects of the quantum vacuum may be taken into ac- 
count introducing a vacuum stress-energy tensor, T™"^, in the Einstein equa- 
tion, which should read 

R,.-^g,.R = k{T,, + Tp'), (1) 
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where R^^ is the Ricci tensor, R the curvature scalar and T^,^ is the stress- 
energy tensor of matter either baryonic or dark plus radiation, k is Sir times 
Newton's constant and we shall use units c = 1 throughout. We assume that 
the vacuum tensor, T^"'^, depends on the space-time curvature, therefore it 
should be a functional of the Riemann tensor, R^y\a (and the metric, g^u-) 
It is plausible to derive the functional from a generalized Einstein-Hilbert 
action 



where F is associated to the vacuum. In general F should be a function of 
the scalars which may be obtained by combining the Riemann tensor, R^uXai 
and its derivatives, with the metric tensor, g^^. Here I shall considers scalars 
more general than the Ricci scalar, whose functions have been extensively 
explored in recent years under the name of f(R)-gravity^,^. 

The action eq.(l2]) may be interpreted as a modification of general rela- 
tivity (RG), that is we may put the tensor T^y'^ on the left side, rather than 
the right side, of Einstein eq.([T]). Indeed it is equivalent in practice whether 
we assume that the quantum vacuum gives rise to some stress-energy to be 
added to the matter one or we assume that GR should be modified by adding 
to i? a term F in the gravitational action. Even if we remain within metric 
theories of gravity, that is we assume that the stress-energy tensor of mat- 
ter produces the curvature of space-time, the arguments leading to the GR 
choice (the standard Einstein equation) are not compulsory, but small mod- 
ifications are compatible with observations. That is we might introduce the 
sum R + F, instead of the Ricci scalar R in the Einstein-Hilbert action. Of 
course there are strong constraints to the form of F from both observational 
evidence and requirements of consistence. 

There are a number of proposals for the gravity of the quantum vacuum 
derived from fundamental arguments involving the quantization of model vac- 
uum fieldspj . They lead to F being a function of scalars like R^^R'^^, R^uXaR'^^'^" 
and OR in addition to R. Most of the proposals have been made for the study 
of cosmology, in particular the attempt to explain the observed accelerated 
expansion of the universe (the "dark energy"), although it is not necessary 
to modify GR in order to explain it0. 

The aim of the present paper is to study the influence of modified GR 
in the structure of Newtonian or weakly relativistic stars. A study of fully 
relativistic stars has been made recently within f(R)-gravity[8j, where the 
authors conclude that neutron stars are not possible (or they require extreme 
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fine tuning) within f(R) theories if these are compatible with observational 
constraints. However the theories studied here are more general than those 
cosidered in Ref.|8]. 

2 Vacuum gravity for weak curvature 

In order to find the most appropriate function F to be put in the action 
eq.(l2]) I shall not attempt to derive it from fundamental arguments, but use 
a plausible phenomenological approach combining arguments of simplicity 
with dimensional considerations. The Ricci scalar has dimensions, L~^, of 
inverse squared length and the theory derived by using it in the action (that 
is GR) is known to give very good agreement with observations for a wide 
range of intensities of the gravitational field (i.e. curvature of space-time.) 
The purpose of this paper is to study the influence of vacuum gravity on 
Newtonian or weakly relativistic stars. Therefore I propose to include in F 
only terms with dimensions not departing too much from L~^. Thus I shall 
assume that F is a sum of terms with dimensions L°, and L~^. For sim- 
plicity I will exclude more complicated terms with dimensions intermediate 
between and L~'^, like L~^logL. Thus I am lead to the following 

F = A + aoR + aiR^ + a2Rf,,R^'' + a3DR + a4V^V''R^, + a5Rf,uXaR'"'^'', (3) 

where means the covariant derivative and □ = g^^V^Viy. The constant 
parameter A has dimensions oq is dimensionless and all remaining coef- 
ficients aj have dimensions L^. There are no other scalar terms with similar 
dimensional dependence so that eq.(|3]) gives the most general F fulfilling 
the requirements of dimensionality and simplicity above stated. However I 
recognize that the stated criterion of simplicity might be questioned. 

In eq.(|3]) the term A will give rise to a "cosmological constant" in Ein- 
stein's equation. It may be relevant in cosmology, but it will have a negligible 
effect in the structure of stars and I ignore it in the following. The term aoR 
may be absorbed in the standard GR term, which amounts to a rescaling of 
Newton constant (the new constant k will be the old one divided by 1 + Oq). 
After this rescaling the gravitational action in eq. contains only the "true 
vacuum polarization", in the words of Zeldovich|9].. The term with as may 
be removed because it gives no contribution to the field equations. The same 
is true for the term with 04 due to the fact that V^{R^u — \g^ivR) = 0, 
which leads to the equality V^V^-R^j, = \^R- Finally the term with 05 may 
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be removed taking into account the well known fact that the Gauss-Bonnet 
invariant, 

g = R^- AR^.R^" + i?^,Aai?^"^", 

does not contribute to the field equations and therefore we may substitute 
4:Rfj,iyR^'^ — R^ for Rfj^i^x^R^^^'^ without changing those equations. Thus I will 
study here the action 



S = ^j d'^xy^ {R + aR" + bR^.R'^") + S^au 



(4) 



which contains just two free parameters, a and h. 

The field equations associated to the action eq. (jl]) may be got from Car- 
roll et al. [To] leading to 



R 



-g^,uR + 2a 



RR^iu — -^QfivR — "^fS^uR + gfiuOR 



+b 



-6 [V,V,Rl + V.V^i?^] = kT^ 



(5) 



where T^i, is the stress-energy tensor of matter. The last term of the left side 
may be transformed taking into account the standard rule for the commuta- 
tion of covariant derivatives of a tensor, that is 



g'-'^V^V^R^x = VuV^Rl + R'^^R 



A per 



^X^fjiva 



(6) 



When this is put in eq.(l5]) its left hand side becomes 



R 



^g^,uR + 2a 



RR 



^g^,,R^ - V^VuR + g^^uUR 



+b 



--^g^iuRxcjR^'' + ^R^lu + -gfiuOR — V^Vj^i? + 2R'^^Rx^ua 



where I have taken into account the equalities 

pcrA p P"^^ R R'^"^ F? F?"^^ R 

and, taking into account that the divergence of Einstein's tensor eq.(IH]) is 
zero. 



V.R 



Act 



The trace of the field equation is specially simple, namely 



(6a + 2b)DR-R = kTj^ = kT. 



(7) 



It is convenient to rewrite the field equation in terms of Einstein's tensor, 
Gf^u, and its trace, G, related to the Ricci tensor, by 



R 



fJ,U 



G 



^dfiuG, R 



-G. 



(8) 



Also I shall write the field equation so that it looks like the standard GR 
eq.([T]), that is 



G 



IcT -I- hT"""-^ 



(9) 



defining 



= k-'{{2a + b)[g^,nG-V^V,G]+a 



GGfj_u — —g^uG'^ 



+b 



Xfii/a 



-gfj.uGxaG^'^ — -g^^G"^ — nG 



}. (10) 



I stress that all results of this paper will be independent of whether we assume 
that T^"^ is a quantum vacuum stress-energy or we consider eqs. ([9]) and f lTOj) 
toghether as a modification of standard GR, maybe with no reference to the 
vacuum. In any case I shall use throughout this paper a language appropriate 
for the former assumption. 



3 Approximate vacuum field in a Newtonian 
star 

Our task is to solve the field equations for a spherically symmetric body, 
with mass M and radius i?o, in a static spacetime with the condition that 
the metric is asymptotically fiat (Minkowskian). The body may be a metallic 
sphere in a laboratory experiment, say like that of Eotvos, the earth or the 
sun, but to be specific I shall speak about a star from now on. I will use 
standard (curvature) coordinates with metric 

ds^ = - exp {13 (r)) dt^ + exp (a (r)) dr^ + r^dVl^. (11) 
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In a static problem of spherical symmetry there are only 3 independent com- 
ponents of Einstein's tensor eq.([n]), namely Gl,Gl,Gg = G^, which are well 
known functions of the metric parameters a and /3[Tl]. Then we have 4 
independent equations including a relation between density and pressure; 
in particular for a star in equilibrium the latter relation is the equation of 
state of matter, i, e, p = p{p). We have also 4 unknown functions, namely 
a (r) , /3 (r) , p(r) and p(r). The solution of these equations is involved in gen- 
eral and here will be solved only for a few particular cases of Newtonian, or 
slightly relativistic, stars. 

Thus I shall consider stars where: 

1" The metric coefficient /3, eq. dlip , is small compared with unity. In this 
case /9 ~ 2$, $ being the Newton potential which fulfils |$| f» kM/Ro << 1. 

2° The matter pressure, p, is small in comparison with the matter density, 
p. Actually this condition is related to the former because we have p/p 
kM/Ro. 

In order to get the 3 components of T^^^{= k~^G^y — T^'/;'^), we shall solve 
eq.dHj). For Newtonian stars the equation may be approximated as follows. 
Firstly we may neglect terms quadratic in G^y. In fact, the terms linear in 
in the right side of eq. lfTOj) are of order akp/Rl, p being the typical 
matter density. In contrast the terms quadratic in G^^ are of order ak'^p'^, 
that is smaller than the former by kpRl ~ kM/Ro << 1- Thus eq.(l9]) may 
be approximated by 



G,y - kT,, = kT;:' ~ (2a + h) [g^^UG - V^V.G] - hUG^,. (12) 

Secondly we may neglect the matter pressure in comparison with the density 
in the interior of of the star, that is we may assume 

T/ = p ^ T, T; = Tl = T; = ^ 0. (13) 

However I shall retain p in some cases for the sake of clarity. Using the metric 
eq.f fTTj) we may write 

where I label 1 the index of the radial coordinate, in order to avoid confusion 
with the coordinate itself. The affine connections Fj^ and Fj^ are of order 
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kM/Ro, whence the terms involving them may be neglected. To the same 
order we may approximate 

g'Xe = 9'^%^ ^ -Ur, 9^' ^ 1, ^ V'G. (14) 

The term DGf^u is more involved, although straightforward, and I shall 
not write it in general. For our case, that is with the metric eq.f lTT|) and the 
approximations eqs.f fT^ , the 3 independent components of the tensor eq.ffT2l) 
become 

bV^Gl + Gl = {2a + b)V^G + kp, 
Ah 9 fir' 

bV^G; + ^{G',-G:) + G; = {2a + b)-^-kp, 

bV'Gl+'^-^^iGl-GD + Gl = i2a + b){^ + ^^)-kp. (15) 

At this moment I point out that, at a difference with general relativity, 
here local isotropy of the matter stresses does not imply isotropy of the spatial 
part of the Ricci tensor. That is, although the Einstein eq.([l]) obviously leads 
to the implication 

Ttr rp0 /~ir /^d , Tir rid 
y> -L Q —r Lt^ Lt^ =7^ iX^ 

this is no longer true for the more involved field eq. (ITOl) (or the approximate 
eq.( IT2|) .) Indeed we may have GJ^ ^ Gq (and therefore anisotropy of the 
vacuum stresses, i.e. (T^°''^Y^ ^ (T^'^'^)g ) even if Tl = Tq. Furthermore 
the anisotropy of the gravitational field, i. e. the inequality G^ ^ Gq, 
might induce anisotropy of the matter stresses, i. e. ^ Tq , but I shall 
not consider that possibility in this paper. Nevertheless the matter local 
anisotropy, ^ , induced by gravitational field anisotropy, G^ ^ Gq, 
might be relevant in strong gravitational fields. 



4 The field outside the star 

The solution of the field equations outside a spherical star, with a general 
f{R), has been studied by several authors [12]. The interest of the problem 
is that it puts constraints on the functions f{R) in order to be compatible 
with known facts in the solar system. Here I shall make a similar calculation 
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for our action eq.(|2]) with the purpose of finding the range of values of the 
parameters a and h compatible with terrestrial and solar system observations. 

Thus our aim is to get the Einstein tensor, C^^, outside the star taking 
into account that = there. Hence the vacuum stress-energy might be 
easily obtained, that is 



P 



In order to solve eqs.f fT5|) I begin obtaining appropriate linear combinations 
of them. If we add the first equation plus the second one plus two times the 
third, we get the trace equation (compare with eq.(I7])) 

G - (6a + 26) VV = kT. (16) 

where I have approximated DT by the flat-space Laplacian of the matter 
density, V^p (see eqs.flT^ and f[T^ ). I shall assume 6a -f- 26 > 0, that is I 
exclude the case 6a + 26 < which, leading to an oscillating function G'(r), 
is unphysical. The limiting case 6a + 26 = will be considered below. 

Eq. (fT6|) in flat space may be solved by Green 's function method as follows. 
In order to simplify the writing I shall sometimes use a dimensionless position 
vector, X, and a dimensionless star radius X defined by 



I I V ^° 6a + 26 

:, z=|x|,A= =, 7 = W — - — , (17) 



V6a + 26 V6a + 26 \j |6| 

where I have introduced also the parameter 7 for latter convenience. Thus 
the fundamental solution for the trace eq. (fT6ll may be written 



VV (x) - / (x) = -47r5=^ (x) ^ / = - exp {-x) , (18) 

X 

Hence the Einstein tensor outside the star may be obtained by integration, 
giving 

G\ + G: + 2GI ^ G{x) = ^[ ^i^exp(-lx-zl) 

47ry|z|<x |x-z| 

47r (6a + 26) r "^""^ \ VG^^E) ' ^ ' 

where 

M* = {Qa + 2hf''^ p{z)4'Kzsin\iz dz, (20) 
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We see that for (small) objects fulfilling Ro « -\/6a + 26, that is X « 1, 
sinh^; ~ z, leading to M* ~ M. This is not the case for bodies such that 
Ro > V6a + 2b. 

The physics behind these result may be better understood if we define a 
new mass parameter 

Ma = M* exp (-X) = (6a + 26) exp (-X) / p{z)Attz sinh z dz, (21) 

Jo 

so that eq.(fT9ll may be rewritten 

A-K (6a + 2h) r """^^ \ V6a + 2h) ' ^ ' 

We see that the vacuum correction eq.(!22l) may be interpreted as if it depends 
on the distance, r — i?o, from the point r where we measure G to the closest 
point in the surface of the star and the correction is due to the mass, Mq, 
contained in some volume of the star most close to the point r. Indeed, the 
factor exp {—X) sinh z in the integral eq.f l2ip effectively restricts the range of 
integration to a region near the surface. In particular for a body of constant 
density, like the earth, or more generally any rocky planet or satellite, eq.f l2T|) 
gives 

M„ 3 , , ,,,, 3 , , ,,,, 3 (6a + 26) 

(23) 

the latter equality being valid for 6a + 26 >> i?Q. In a star, like the sun, the 
ratio Ma/M would be much smaller because the surface density is smaller 
than on earth and the radius much bigger. 

In order to proceed with the calculation of the Einstein tensor outside 
the star I subtract from twice the first eq. IlT^ the second one and twice the 
third. This gives 

6V' {2G\ -Gl- 2G'l) + 2G\ - G', - 2Gl = 2p, (24) 

again neglecting p << p. The solution outside the star is (compare with 
eq.(l22])) 

2Gl - G: - 2Gl = exp ^-'-^] , (25) 
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where 

Mb = {-hf^ exp (-7X) / p{z) sinh zAnzdz. (26) 

Jo 

I have assumed 6 < because a positive h would lead to an unphysical 
oscillating function. For a body with constant density this leads to 

t = ^ 11 + exp (-2,X)\ - ^ II - exp (-2,X)] - H « L 

(27) 

the latter equality valid for |6| « R^. 

From eqs.f l22|) and f l25l) we get outside the star 

pvac ^ j^~iQt ^ ^_ i^y {X-x)+ 2Mb-f^x^ exp (7X - jx)] , (28) 

pr + 2pr = -A:-' (G; + 2G^) 

= ^ [M67^x^ exp (7X - 7a;) - MaX^ exp (-x)l . (29) 

In order to obtain separately the two different pressures, p'^"''^ and pg"''^, I 
proceed as follows. The vacuum stress-energy tensor eq.f llOp is divergence- 
free as may be easily checked. Actually this property is a consequence of 
deriving the field equations from an action functional. As the Einstein tensor 
has also zero divergence, eq.(|9]) shows that both the matter and the vacuum 
stress-energy tensors are separately divergence-free. The divergence of the 
vacuum tensor gives the relation 

J^vac vac ^^vac 

dr r 

where /3' is the radial derivative of the metric coefficient /3. The second equal- 
ity follows from the fact that, in Newtonian stars, is of order kM/Ro << 
1, whilst p^"''^, p^"''^ and pg"''^ have the same order. Hence we obtain 

i. (^3p-e) ^ ^2 ^pvac ^ ^^vac^ ^ 

which, taking eq.( 12U]) into account, gives 

= ^ [Ma {x + 1) exp {X-x)- Mb (7X + 1) exp (7X - 70;)] , (31) 
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pT" = [^b (7^^^ + 7^ + 1) exp (7X - ^x) -Ma{x^ + x + 1) exp (X - x)] . 

(32) 

where the labels eqs.( IT7j) have been used. An integration constant has been 
fixed so that r^p^"'^ — )• for r — )■ cxd in order that the absolute value of Pr 
never surpases p, eg. (1251) . 
The limiting case 

3a + 6 = 0,-6 > 

cannot be studied by the procedure leading to eqs.(!28l) to (136!) but its trace 
eq.( ll6p is rather simple, namely 

G = kT =^ T'""''^ = 0, in particular G = outside the star. 

Also eqs.lTT^ are simple and we obtain outside the star (here Ma = M^) 



„vac 



Pr 



ISirar \ a/3^ 
-k G,, = -- — - —= + 1 exp 



67rH VV3a / V v3a 



10 Ma f r \ f r-Ro 

-k ^Gq = — — ^ — H 1= + 1 exp 



127rr3 \3a JSa J \ V^a 



We see that these equations are the limit of eqs. ( 128|) to f p2|) when b — )■ —3a < 
0. It is interesting that, in this case, the vacuum field looks like a radiation 
field because its stress-energy tensor is traceless. However that radiation field 
is not isotropic, that is p^"""^ 7^ Pe"''^, Pr"''^ being negative and pg"^ positive. 

Also the case 6 = cannot be studied by the procedure leading to eqs.( l28|) 
to (136|1 but eqs.( IT5|) may be easily solved taking the trace eq.( !22|) into account. 
We get 



, vac r-e 2a rfG kMa ( r\ r \ ( r - R^ 
kpf = Gg = -G — = — - — - — + —= + 1 exp ' 



3 r dr \Qa V6a / V VOa 

(35) 
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It is interesting that, in this case, the mean vacuum pressure is negative 
fulfiUing p2L^ = -Ip^--. 

Finally in the particular case 

2a + 6 = 0, 

eqs.( ]2S]) to (15^ hold true, but it is interesting to write them explicitly because 
they are specially simple, that is 

G = G\ = kp^'^'^ = ^ exp (-'-^) , G: = pr = Gl = pT = 0. 

In summary the vacuum density and stresses outside the star, derived 
from the action eqs.(Il]) , are given by eq. (l28l) and (!29|) . or the appropriate 
limits, for any values of the parameters a and h fulfilling 

3a > -6 > 0. (36) 

For values violating these inequalities the solution of the field equation out- 
side a Newtonian spherical body would give an unphysical oscillatory be- 
haviour. Of course the values a = 6 = correspond to general relativity 
(without vacuum field.) 

Now I shall calculate the coefficients a (r) and /3 (r) of the metric eq.f llip 
taking into account the well known relations 

exp (-«) = 1 - - / u'Gl iu) du, 4- = P' = re'^G:. (37) 

r Jq dr r 

From the former it is straightforward to obtain the function a (r) outside the 
star. We get, taking eq.( l28|) into account, 

exp(-«) = l- — , 38 

Anr 

where 

6a (x) = i^[(x + l)exp(X-x)-(X + l)] 

+ [{ix + 1) exp (7X - 7x) - (7X + 1)] . (39) 

We see that the dimensionless function 6a represents the correction to the 
GR (Schwarzschild) exterior solution of a spherical star. 
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In order to get the function /3' (r) I start expanding the metric coefficient 
exp a in powers of the gravitational constant fc, a parameter which may be 
considered small because kM/Ro « 1. That is 

exp a = 1 + - — {l-5a) + — ^ 1 - 5a + O (k^) . 

Inserting this in the second eq.( l37|) we obtain as an expansion in powers 
of k, that is 

(5' = -— {l-5a + 



47rr2 \^ M 



1 - 5a) (^1 - 5a + -^-^ + O (A;^) , (40) 



where 5a was given in eq.f l39|) and 

1-Sa + ^^ = ^[{x + l)exp{X~x) + {X + l)] 

+ ^ [-4 iix + 1) exp (7X - 7x) + 2 (7X + 1)] . 

This leads to 

+^ l(7A' + 1) - 2 (71 + 1) exp (7A' - + O (i^) . (41) 

Getting the term proportional to k^ is straighforward, but the resulting ex- 
pression is involved and I will not write it explicitly. 

Hence the parameter /3 may be easily obtained by means of an integration 
of eq.f HOj) with the condition /3 (r) — ?■ when r — )■ oo. The result may be 
written 

/'=-i^(l + W-|^(l + W + 0(t3). (42) 

where the function 5/3^ (r) is the correction to the GR (Schwarzschild) solu- 
tion to lowest order in k. It is 

5/3i = ^ [(X + 1) + exp(X - ^)] + ^ [2 (7X + 1) - 4exp(7X - 7^)] . 

(43) 
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The term 5/32 is involved and I will not write it explicitly.. 
Eq.f l42p to order 0{k) is, therefore, 
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kM 



3M 



3M 



kM, 



127rr 



— exp(X — x) H exp(7X — 7a;) 



Svrr 



(44) 



As is well known the term linear in the gravitational constant, k, in the 
expansion of the metric parameter f3 equals twice the Newtonian potential. 
Therefore eq. f l44p shows that the theory resting upon the action eq.(|l]) pre- 
dicts: 1) a correction to the mass appearing in the Newtonian potential (the 
term within square bracket), and 2) two non- Newtonian potentials of Yukawa 
type, one of them attractive and the other one repulsive. 



5 Constraints on the parameters a and b 

It is known that in f(R)-gravity the coefficient of the term i?^ (noted a in 
this paper) cannot be greater than a few square milimeters in order not to 
contradict laboratory experiments [4j. However this is not the case in the 
generalized theory here studied. In fact, let us compare the Newtonian field 
intensity, g = kM/ (87rr^) with the field intensity eq. pT]) predicted by our 
theory, which corresponds to half the term linear in k of eq.f l4ip and I shall 
label g + 6g. We get 

- = |§[(x+l)exp(X-x) + (X + l)] 
g SM 

+^ [(7^ + 1) - 2 (70; + 1) exp (7X - 7x)] . (45) 

In the most interesting case of a body with constant density we may put 
eqs.(l23!) and (!27|l in eq. (H5|) and we get 

6g X coshX — sinhX ^, , , ^, , 

— = 77^ (x + 1) exp (-x) + (X + 1) exp (-X) 

g X-^ 
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This expression is small for both small and large X. Indeed for X ^ x << 1 
we have, expanding in powers, 

(47) 

For X >> 1 and x — X « X (i.e. r — Ro « Ro) we obtain 

5g 1 + 27"^ X .^^,2x . , 
— - + ^exp(X-x)---^exp(7X-7x) 

VQa + 2b + 2J\b\ V6a + 2b f r - Ro \ ( r - R^ 
H exp -^=^= ^ — exp 



Ro Ro V v6a + 26/ Ro \ a/|&| 

It is a fact that precise measurements have been made only for the grav- 
itational field of objects with sizes of meters or smaller (e.g. the Eotvos 
experiment) or for celestial bodies, in particular the earth[T3]. Thus if both 
y/a and have values between hundred meters and a few kilometers both 
eqs.f H7|) and (jUj) predict violations of Newtonian gravity of order 10~^ or 
less. The two parameters should be large enough because if |6| is small then 
the ratio 5g/g becomes large in eq. fH^ and the theory here developed is 
refuted by laboratory experiment s|13j. It is worth to remember here that 
during the late eigthies of the XX Century there were some claims about the 
existence of non-Newtonian gravity (the "fifth force" ) , although a reanalysis 
of the experiments has lead to a consensus that several uncertainties were 
not taken into account and the experiments are actually compatible with 
Newtonian gravity within errors [H]. In any case I should mention that some 
of the experiments apparently showed the existence of two non-Newtonian 
potentials of Yukawa type, one attractive and the other one repulsive, as in 
our eg. (1481) . For instance in a tower experiment Eckhardt et al.p2] reported 
a non-Newtonian gravity in the form 

r 

^at J \ ^rep , 

with values 

Erep - Sat ^ 0.007, E^t > 0.03, Xrep ~ A,^ ~ 100 m. (49) 

These values violate our eq.(H8ll which predicts 

A 



1 + e^t exp ( -— - Erep exp 

Ant / \ \ 
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In constrast this equality is compatible with known bounds provided that 
i/a < 50 km and e < 0.005 (see Fig.l of Ref.|14j.) I conclude that tests of 
our theory in the earth surface would require experiments with errors several 
orders smaller than those typical of "fifth force" experiments. 

In summary there are two ranges where the parameters are compatible 
with all performed experiments on Newtonian gravity, namely either both 
a/o and are less than a few millimeters or both have values between 
hundred meters and a few kilometers. But I shall point out that, if y/a and 
^/^\ are in the latter range, the correction eq.( H6|) may be quite important 
for mass concentrations with sizes of a few kilometers, like mountains. This 
may lead to experimental tests of the theory, but I shall not discuss them in 
this paper. 

Even if the discrepancies between Newtonian gravity and the predictions 
of eq. poj) are too small to be detected, it is interesting to see whether the 
theory here developed predicts corrections to Newtonian gravity in measur- 
able violation of general relativity. The standard comparison is made using 
an isotropic metric rather than eq. ( 1TT|) . However it is possible to make the 
comparison also with the latter metric [TT]. For a spherical body the coeffi- 
cients of the metric eq. (fTT]) may be expanded in powers of the gravitational 
constant, k, in the form 

kM ^ kM ^ ^ k'^M^ , , 

expa = 1 + r^— + exp/3 = 1 - A— + (C - + - (50) 



General relativity predicts \ = = i] = \. Furthermore, recent measurement s [18 
give the bound ry/A = 1 + (2.1 ± 2.3) x 10"^ 

In order to make the comparison with the GR prediction eq.f l50p we need 
the expansion 

exp/3 = 1 + /3 + + O (A;3) 

kM k'^M'^ 
- 1 - ^(1 + + (1 + 25/3, + {5^r - 5/3,) . (51) 

Now the comparison of eq.flSOl) with eqs.(l38l) and ( ISTI) gives (to lowest order 
in k) 

X = l-6a,7] = 1 + (5/3i, C = 1 + 35/3i + {6/3^^ - 6/3^ -6a- 5a5P^. 

With arguments similar to those used in the analysis of eq.( l45l) we get A = 
7] = C = 1 with errors smaller than 10~^ for the Earth gravity and smaller 
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than 10~^ for the gravity of the sun, provided that the parameters a/o and 
are less than a few kilometers. The reason for the difference between 
terrestrial and solar gravity lies in that, according to our predictions, the 
quantities 5a, 6(3i and 6(32 ^^re proportional to the density near the surface 
(much smaller in the sun ) and inversely proportional to the radius (100 times 
greater in the sun). In particular I get 



1- 



2M, 



3M 



{'jx — 1) exp (7X — 7x) + -^TTjix + 2) exp {X — x 



< 2.10"^ 
(52) 

in agreement with observations [18] . 

For the gravitational interactions of bodies with sizes of a few meters or 
smaller, used in laboratory experiments, the ratio rj/X is of order R^/a < 10~^ 
if ^/a and ^/\b\ are larger than about one kilometer. 



6 Hydrostatic equilibrium 

The structure of a spherical star in equilibrium may be obtained from eqs.(l9]) 
and (fTOl) plus the equation of state. The solution of these equations is involved 
and furthermore the fact that they are fourth-order shows that we need four 
(initial or boundary) conditions, which may give rise to some ambiguity. I 
will not attempt solving exactly those equations in this paper, but I shall 
sketch the solution of the approximate eqs.( IT5ll inside Newtonian or slightly 
relativistic stars. 

I shall begin solving the trace eq.f fT6|) , via eq.f fTOj) , in the interior of the 
star, that is x < X. I get 

k 

G (x) = — / p (z) [exp {— \x — z\) — exp {—z — x)] zdz. (53) 
2a; Jo 

For a star where X >> 1 we may make the following 3 approximations: 1) 
neglecting the second exponential in comparison with the first one (except 
near the origin), 2) approximating p{z) by 

{z - xf , 

and 3) extending the z-integration to the interval (—00, 00) . Then the integral 



p{z) ^ p (x) + 



dp 

dz 



, 1 



(Pp 
d? 
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is trivial and we get 



GW = ..(T + T"n-*^(.W + |^ + |g). (54) 

This same result is obtained near the center although the approximations 
involved should be different. It may be realized that the last two terms of 
eg. (15^ are of order = (6a + 2b) / « 1 and the terms neglected are 
of the order of X~^. 

The same result may be obtained if we combine eqs. ( !T6|) and (Q , the lat- 
ter with substituted for □, and work to the same order of approximation. 
In fact we get 

T"""^ = (6a + 26) (T + T'"^") ~ (6a + 2b) V^T ~ (6a + 2b) V^p. (55) 

The agreement between both calculations, leading to eqs. ( 15^ and (!55|) , re- 
spectively, reinforces the validity of solving eqs.( IT^ by approximating by 
A;T^ in all terms which are linear in the small parameters a or b. I point out 
that this approximation was not valid in the study, made in Section 4, of the 
exterior of the star because = there. Thus eqs.f fTSj) lead, in the star 



interior, to 



- {Tnl^2a(-^ + f^)^2a(^-p' + p" 
\r ar dv^ J \r 

pvac ^ _(^J^vacy ^_^2a + b)-p', 

r 

pr ^ -(T-^)^^-(2a + 6)(^V + p")- (56) 

These results are consistent with eq.( l55|) for the trace T'"^'^. It is remarkable 
that the vacuum density in the interior of the star depends only on the pa- 
rameter a, and not on b, within our approximations (in particular neglecting 
the matter pressure in comparison with the matter density). If —6 < 2a the 
mean vacuum pressure is negative in the interior of the star. 

I shall point out that eqs. ( 15^ include all vacuum effects to our order 
of approximation. In particular we should not add to eq.dSSD the effect of 
the Yukawa-type contribution to the Newtonian potential (see eqs. pUj) ). It 
is remarkable that the density of the vacuum contribution is negative in 
the central part of the star because both dp/dr and d'^p/dr'^ are negative 
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there. However the total density is positive everywhere when eg. f l56l) is vahd 
because we have assumed that \p'"°''^\ << p in deriving it. The vacuum density 
is positive near the surface of the star and the vacuum contribution to the 
total mass is zero, that is 

f 2a (-^ + ^] Anr^dr = Snar^"^ < 0, M"'^'^ = m™'^ (Ro) = 0, 
Jq \r dr dr^ J dr 

(57) 

provided that dp/dr is zero at the star surface. In particular this is the 
case in all polytropes with 7 < 2 (see next Section.) I stress that all these 
properties are valid only for large spherical static bodies, where eqs.( l56l) are 
a good enough approximation. I shall point out that, strictly speaking, the 
total mass of the star should include the mass of the vacuum contribution 
outside the star (with density given by eg. (1281) ) so that M'""''^ as defined in 
eq.f lFr|) is only the mass in the interior of the star. The external mass due to 
the vacuum may be obtained by integration of the density eg. (1281) giving 

the latter eguality being valid for a celestial body of constant density like the 
earth. 

The three independent components of the Einstein- type eg. ([9]) become 

Gl = kp + 2ka ( -p' + p") , = -kp + (2a + b) k-p\ 
\r J r 

Gl = -kp+{2a + h)k(-p' + p") , (59) 



r 



where p' = dp/dr, p" = d'^p/dr^. Egs.f l59|) will be the starting point for all 
calculations of structure of stars to be made in the following. 

The condition that Einstein's tensor GJ^ is divergence-free leads to the 
following hydrostatic eguilibrium eguation[T9] 



dr r \ J 2 Anr"^ — krm'^ff 

where I define 

p-ff = p + p--, p^ff =p + p--, pl^f =p + pf 



m 



= An r'^p^^^dr. (61) 
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Putting eqs.dSS]) into eq. we get 



_ k,n + 4,rV - 8. (a + b) r't (^^^_^,\i±^ 2a^\ . (62, 



dr 2 47rr2 - krm - Snakr^^ \ r dr dr'^ 

I point out again that this equihbrium equation might be seen as either a 
modification (to first order in a and b) of the general relativistic equilibrium 
equation, due to the change from the Einstein-Hilbert action to the action 
eq.(|2]), or as an effect of the quantum vacuum. Both interpretations lead to 
the same physical consequences. 

In going from eq. ( 160|1 to eq. (!62|) it is interesting the cancellation 



which actually derives from eq.( l30|) . A fortunate consequence of this cance- 
lation is that eq. (l^2]) is an ordinary differential equation of third order (in 
the variable m(r) ), rather than fourth order, when p is written in terms of 
p, e. g. via the equation of state. Then we may use the initial conditions 
m (0) = 0, [dp/dr]^^Q = and only one more condition is needed, e. g. 
fixing the central density, p (0) . Thus the solutions of eq.f l62|) consist of a 
one-parameter family as in the standard theory of stars in equilibrium. 

Assuming that both the relativistic and the vacuum corrrections are 
small, eq.( l62|) may be approximated as 

dp km / ^„ . ,s 

— = -p [1 + GR + vacuum) , (64) 

dr Snr^ 

47rr^» km p 

GR = - + - — + -, 

m Airr p 

, ^.%'Kr'^dp 1 dp ^ Id^p 

vacuum = — (a + o) ; 26 — h 2a-— -r. 

m dr rp dr p dr^ 

It may be realized that the first factor on the right side represents the Newto- 
nian approximation, the 3 terms labelled GR are the well known first order 
corrections of general relativity and those labelled vacuum are corrections 
due to the vacuum stress-energy. The former (latter) are of order kM/R 
(order a/R^) with respect to the Newtonian approximation. I have ignored 
terms which are both vacuum and GR corrections, that is of order akM/R? 
with respect to the Newtonian approximation. As may be seen, and it is 
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well known, the three GR corrections are positive, that is every one con- 
tributes to the increase of gravitational effects. In contrast the second term 
of vacuum is negative (remember that 5 < and dp/dr < 0) and the third 
one is negative in the central region of the star because d'^p/dr'^ is negative 
there. The first vacuum term would be positive (negative) everywhere if 
|6| < a {\b\ > a) . Finally the third term is positive near the surface. Conse- 
quently from eq. flM|) no conclusion seems possible about whether the vacuum 
corrections increase or decrease the effects of gravity in comparison with the 
Newtonian approximation. 

There is, however, a simple argument which gives the answer. The grav- 
itational interaction energy, in the Newtonian approximation but with the 
vacuum correction included, is given by 



where we have taken eq. ( 157|) into account and neglected the vacuum density 
in the exterior of the star (but see below). Also I have performed an inte- 
gration by parts in the second equality assuming that p = at the star's 
surface. The reason why the vacuum correction is positive derives from the 
fact that m^"'^(r) < (see eq.f l57|) ) so that the vacuum contribution makes 
the Newtonian gravitational interacction energy less negative. It is as if there 
were a short range repulsion given by 



Actually the vacuum density outside the star gives another contribution to 
the gravitational interaction which is negative. In fact, taking eq.( !28|) into 
account we have 



the latter equality being valid for a celestial body of constant density like 
the earth. Considering the two terms of eq.f lB^ plus the terms eqs.f lSS]) and 



V 




3kaS^ (ri - rg) . 



rvac 
ext 




(66) 
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(IMj) we conclude the following. For the earth the Newtonian gravitational 
energy, V, and the vacuum term, M'^^^, may be of similar size (their ratio 
is of order kM/ s/a) , whilst the other two terms are much smaller. However 
in star, where the surface density is very small, the dominant term may be 
V , followed by V'""''^, whilst the other two terms are much smaller than the 
latter (for the proof in the case of white dwarfs see next Section). 

The relative importance of the vacuum and GR corrections may be esti- 
mated from the ratio of the first vacuum correction to the first GR correction 
in eq. (!M|) . It is 

vacuum 2a + bdp a p , , 

~ ■ ~ — X -. (67) 

GR rp dr R^ p 

In Newtonian stars the first factor on the right side is small, as was shown 
in the second Section, and the second factor is large. Thus no conclusion 
may be reached (in any case both GR and vacuum corrections are small). In 
compact objects like neutron stars the second factor is of order unity and the 
first one might be large, which suggests that the vacuum corrections could be 
very relevant. However in these stars many of the approximations leading to 
eq.f lM|) are not valid and a calculation starting from eqs.dH]) and ( ITU]) would 
be necessary. 

As an illustrative example in the following I study the effect of the vac- 
uum corrections on the internal properties of the sun. I shall consider, in 
particular, the change in central pressure due to the correction terms given 
in eq. fl64p . The change in the central pressure will be 6p (0) where 

op [r) = — ——dr = — / -Tipdr \vacuum\ (68) 
J, dr 87r J,. r^ 

/"^^ Z"^® m dp^ m d^p^ 

= —ia + bjk p—dr — bk / - — -—dr + ak / - — -—^dr. 

J J. dr J J. Airr-^ dr 47rr^ dr'^ 

An integration by parts of the third term gives an integrated part, a term 
which combines with the first one and can be integrated easily, and a new 
term which combines with the second. Thus we have (assuming that p' {Ro) = 
0) 

6p (0) = (a + ^) A: [P (0)]' + (2a - b) k -^^dr. (69) 

The integral in the second term may be performed numerically using data 
of the calculations made by Bahcall et al. [20], but for our purposes it is 
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enough to find an upper bound. Taking into acount that m/r^ < ^np (0) and 
dp/dr < 0, the value of the second integral lies in the interval [—^ [p (0)]^ , O) . 
Thus we get 

la + ^b^ k [p iO)f < 6p (0) < (a + (0)]' , 
which shows that 

\Sp{0)\ ^ (10 km)^ X 7.4 x lO^^^cm/g x (l56g/cm^)^ = 1.8 x 10"^^ g/cm^ 

where I have estimated the parameters a and |6| in about 10 km. We see that 
6p{0) may be negative (e.g. for |6| > 2a) and positive (e.g. for \b\ < 2a/5.) 
The central pressure is (in units of mass/volume) 

p = 1.19^^ = 2.6 X 10-^ g/cm^ 

where fcs is Boltzmann's constant, itlh the mass of the hydrogen atom and 
B the the temperature. The factor 1.19 is the mean number of particles per 
baryon, which derives from the chemical composition in center of the sun, 
that is 35% hydrogen and 65% helium. Hence we get 

6pi0) _6QiO) 



p{0) G(0) 



10-^ (70) 



We see that the effect of the vacuum stress-energy in the structure of the sun 
is negligible. In particular it gives no change in the prediction of neutrino 
emission from the center of the sun. 



7 Stability of 7 = 4/3 polytropes. Application 
to white dwarfs 

In the following I shall apply the theory resting upon the Einstein-type 
eqs.dSnD to the study of equilibrium and stability of 7 = 4/3 polytropes, 
that is Newtonian stars with equation of state of the form p = Kp"' , 7 = 4/3. 
In order to make more transparent the comparison with the literature, in 
this section I shall use Newton's constant, G, rather than k = SnG. No con- 
fussion should arise with the Einstein tensor, which will not be mentioned in 
this section. 
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It is well known that polytropes are stable if 7 > 4/3 and unstable if 
7<4/3. If 7 = 4/3 the stability depends on small corrections which there- 
fore become relevant. In particular, general relativistic corrections produce 
unstability, although other corrections may compensate for that. Stars which 
may be treated as 7 = 4/3 polytropes are white dwarfs and supermassive 
stars. Both are stable when the internal energy is large enough, but become 
unstable after radiating a fraction of that energyPTj. In both cases it is as- 
sumed that a source of unstability are the relativistic corrections, although 
in some white dwarfs also neutronization of the core may induce unstabilil- 
ity. Here I shall study white dwarfs but not supermassive stars. Indeed the 
latter are hypothetical and there are no observations on them. In addition 
the possible corrections due to the quantum vacuum, as they are proposed 
in this paper, are probably small because typical dimensions of supermassive 
stars are far larger than the parameters a/o and -\/— 6. 

White dwarfs are formed from ordinary stars after a period of cooling. 
The theory here considered corresponds to stars sufficiently cold to be treated 
as if the temperature is zero Kelvin. Also I will consider only massive white 
dwarfs because small ones are well approximated by 7 = 5/3 polytropes, 
they are always stable and have no interest here. The theory of white dwars 
in our approach might consists of solving the hydrostatic equilibrium eq.( l62|) 
with the equation of state 

, , dp du 

p = niHn + u{n) =^ p = n- p = n— , (71) 

dn dn 

n being the baryon density, p the mass density and p the pressure. The 
mass mu is close to that of the hydrogen atom (its precise value depends on 
the chemical composition of the star.) The function u (n) corresponds to a 
7 = 4/3 polytriope with small corrections, that is 

u{n) = mimunY'^ + 5u{n) ^ p ~ Kp^'^, u - 3Kp'^/^ (72) 

where is a constant. In the Newtonian approximation eqs. fl62l) and flTTl) 
become 

dp Gm dp 3 Gm{r) 2/3 /^qn 

If we introduce in eq.f TTSj) the appropriate corrections to lowest order we 

get 

dp 3 Gmir) . ,0/3/1 ^r-. \ 

-i- = — ^pirV^-^ (1 + GR + vacuum + eos) , (74) 
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where GR and vacuum were given in eg. f l64p . The additional term, eos, 
corresponds to the modification, 5u{n), of the polytropic equation of state 
mentioned in eq. fl72|) , which I will not write explicitly here. Only the central 
density, p(0) = p^, is needed in order to fix one solution of eq. (l7il) (see 
comment after eq. fl63|) ). Thus for every value of the solution of eq.flT^ 
provides the functions p{r), p{r) and n(r). Hence we could calculate the 
baryon number, N, and the mass of the star, M^^^ , as functions of p^, using 
the well known expressions 



Jo a/1 — 2G'm^-/'-^ (r) /r Jo 

(75) 

Ro being the radius of the star and p*^-^-^ the sum of the matter and vacuum 
densities (see eq.( !6T|) .) The upper limit of the integral giving the total mass, 
M^^^ , is oo, rather than Ro, because we should include the mass associated 
to the vacuum density both inside and in the neighbourhood of the star (see 
eq.(l58l) .) Actually the quantity of interest is the binding energy defined by 

E = M^^f -tuhN. (76) 

which for Newtonian stars (i.e. without the corrections GR, vacuum and 
eos) becomes 

E = 3K p^/^4nr^dr-G / mpAnrdr. (77) 
Jo Jo 

In summary eqs. f l74p and f l76p provide a one-parameter family of solutions of 
equilibrium, and the question is which of such solutions correspond to stable 
equilibrium. 

A standard method to study stability is to start from an equilibrium 
configuration and perform the transformation r — )■ Ar, p — i- X^^p which leaves 
the mass unchanged, that is M — )• M. The binding energy, E, becomes 
a function of A and the configuration is stable if the function E (A) has a 
minimum for A = 1. Thus the equilibrium and stability of white dwarfs may 
be studied as follows. The solution of the Newtonian eq. fl73|) for a central 
density p^ provides the density p (p^; r) as a function of the radial coordinate. 
We may assume that this function is a good approximation for the solution 
of eq.f in|) with the said central density p^. Thus we may use it in order to get 
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the baryon number, N, and the binding energy, E, via eqs. (175]) and (1751) , as 
functions of the central density. Now we consider star configurations out of 
equilibrium, but close to the one given by the function p (p^; r) , by performing 
the A transformation above stated. Thus we get the binding energy, E (p^. A), 
and the baryon number, (p^. A), as functions of p^ and A. They correspond 
to either equilibrium configurations (i.e. fulfilling eq.( 175]) ). when A = 1, or 
non-equilibrium configurations, when A 7^ 1. Now rather than taking p^ and 
A as independent variables we may use and A, finding the binding energy 
as a function E{N, A). It is common to change the variables to M = muN 
and p* = ^^^Pc) so that the function should be written E = E{M,pl), or 
simply E = E{M,pJ. Then the conditions of equilibrium and stability are 

dE{M,p^)/dp, = 0, d'EiM,p^)/dpl > 0. 

But we see from these arguments that we must guarantee that M is indeed the 
number of baryons N times a constant tuh- The point is important specially 
for the calculation of the GR corrections [21]. 

Both eqs.( !75|) and (!76|) contain three corrections to Newtonian theory. 
The first one derives from general relativity and for white dwarfs is of order 
GM/Rq ^ 10~^. The second is the correction for the vacuum contribution 
developed in this article, which is of order a/R^ 10~^ (remember that 

~ a < 10 km and i?o ~ 4 x 10^ km.) The third one is due to the deviation of 
the equation of state from the polytropic one. The 3 corrections are small (the 
latter for massive enough stars), so that we may calculate each of them as if it 
was alone, and add the corrections at the end. Calculating the GR correction 
is delicate due to the fact that the proper volume element is not Anr'^dr. This 
problem however does not appear in the vacuum correction which may be 
treated, in our approximation, as if we solve a purely Newtonian problem 
with a density and a pressure modified by the vacuum contributions as given 
in eq.( l56p . I will take the corrections GR and eos from the literature [21], 
that is 

E = B [KfJl - Mpl'^ + CMp-^'^ - DM^'^pI/^ (78) 

where B,C,D and the Chandrasekhar mass, M^^, are positive constants (I 
use the notation of Ref.[2Tj.) 

In order to calculate the vacuum correction I begin writing the binding 
energy, eq.( 176|) , taking M*^'^-'^ and from eq.( 175|) . The result should be writ- 
ten to first order in G and K (which excludes general relativistic corrections. 
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which are of order or GK). The quantities p^^^ and m^^^ are taken from 
eqs.( 1SB]) and ( IFT]) respectively, which includes the vacuum correction to first 
order in the parameters a and h. This gives 

r-Ro 

E = Arrr'^dr [(p — rrifjn) — Gmr^^mHn\ (79) 

Jo 

"Ro 



2 dp , d'p\ ^^^_dp^ 



2a I -- — h — — — SirGar—mHn 
r dr dr^ I dr 



oo 



+ / Pext4vrr dr, 
Ro 



the first term of each one of the first two integrals being the internal energy 
and the second one the gravitational energy. The vacuum correction is given 
by the last two integrals in eq.f l79p . I shall start estimating the third integral. 
The density outside the star is given by eq. (|25l) with Ma and M^, as in eqs.(l^ 
and fl26l) , respectively. It is known from the Lane-Endem solution of the 
Newtonian eq.f l73l) that the density near the star surface is of the form 

/^N 3/2 

p{r)^Cpl{-^ {Ro-rf, 



where C is a numerical constant. If this is put in the expression of Ma, 
eq.f l2T]) , we get 

/ \ 3/2 

Ma ~ 247tGpI - (6a + 2bf^ Ro. 



Similarly we obtain 



K 



\ 3/2 

M,^2AnCpi[-\ i-bf'Ro. 



If we put these expressions in the vacuum density outside the star, eq.f l28|) , 
we get 

^ p^^A^r'dr = SnGpl j^-j [(6a + 2bf + 2 {-bf] R^ 

The relevant result is that the integral is of order a^, therefore negliglible 
in comparison with the second integral of eq.f l79|) which gives therefore the 
main contribution to the vacuum correction E""""^ (the first integral in eq. f l79|) 
is the Newtonian binding energy). 
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The function p (r) is the density of a 7 = 4/3 polytrope and, for the 
calculation of the integral, we may approximate m//n(r) ~ p (r) . We get 



— GpAnr — 

dr"^ dr 



(80) 



The first term is zero (it equals M''"^, see eg. (1571) ) and the second one leads, 
after an integration by parts, to (see eq. fl65l) ) 

E"""" ~ UiraG / pHnr^dr. (81) 
Jo 

It is interesting that the vacuum correction depends only on the parameter 
a, but not on b. 

The integral eq. flSTl) may be performed numerically in terms of the central 
density, p^ , and the Lane-Emden variables ^ and 9. We get 

E'""' = 12naGMp^^f\9'{^^)\ O^^^d^ 23 AaGMp^. (82) 

Jo 

Adding the vacuum correction, eq. (1821) , to the standard expression for the 
energy we obtain from eq. flTOl) . 



E = B (^M^/f - M2/3 j Mpy^ + CMp^^/^ - DM^/^pI/^ + FMp^, (83) 
with 

F = 23.4aG. 

The first term in eq. fl83p is the Newtonian energy, the second is due to the 
departure of the equation of state from a 7 = 4/3 polytrope, the third one 
is the correction of general relativity and the last term is the vacuum correc- 
tion. For a given mass, the central density of equilibrium is obtained when 
dE/dp^ = 0, which leads to 

B {m^I - j _ Cp;2/3 _ 2DM^/Vy' + 3Fp^/3 = 0. (84) 

Taking into account that the parameters 5, C, D are positive, it is easy to 
see that if F = no value (positive) of p^ fulfils eq. flH^ whenever M > Mch, 
that is equilibrium is not possible. In sharp contrast a value of p^ fulfilling 
eq.f lM|) exists for any mass if F > 0. In fact the quantity in the left side of 
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eq.( 151|) approaches +00 when — )• 00 and —00 when — )■ 0, so that it is 
zero by continuity for some positive value of p^. 

For a mass, M, well below the Chandrasekhar limit the terms with D and 
F may be considered small in comparison with those with B and C, whence 
the central density of equilibrium is 

Pc^Po + —^Po - i^pT. Po = y-^j - M^") (85) 

When M approaches the Chandrasekhar limit, Mchi eq.f lH^ diverges and the 
approximations leading to it are not valid. I shall not attempt to solve eq. f l84p 
for M close to Mch-, which would be cumbersome and not very interesting. 
When M » Mch an approximate solution of eq.( l8^ is again possible, 
leading to the simple expression 

3/2 

^ 8.7 X (^m2/3 _ mHI) . (86) 

However the central density, p^, obtained this way is so big that the approxi- 
mations leading to it do not apply (e.g. for M / Mch = 1-01 we get p^ ^ 10^^ 
kg/m^, greater than neutron star densities.) 

Stability requires that (fE/dpl > 0. When dE/dp^ = the stability 
condition leads to 

Cp-^ - DM*/^ + 3Fpy3 > 0. (87) 

If F = 0, equilibrium is stable for any mass below a limit smaller than but 
close to the Chandrasekhar mass. If F > stars may be in equilibrium for 
any mass, as said above. That hydrostatic equilibrium is always stable, i. 
e. for any mass, M, and for any value of F > 0. In fact eq. flH^ may be 
rewritten 

3FMp2/3 = B (^M2/3 _ M^'Jf) M+CMp-2/3+2F)M^/Vc^^ > DM^'^-Cp-\ 

so that eq. flHTj) holds true (taking into account that i?,C, D > 0, M > 
Mch-) Actually the derived results for stars above the Chandrasekhar limit 
are rather academic because for the probable value of F, that is when the 
parameter ^/a lies below the limit derived in Section 4, central density of 
equilibrium of stars with M close to or larger than Mch would be so large that 



B [hPl^ - Mill 
3F 



29 



the star becomes unstable against neutronization before reaching hydrostatic 
equihbrium. 

In summary, vacuum corrections following from the action eq.(|l]) might 
give rise to dramatic changes in the equilibrium and stability of white dwarfs, 
but these changes would produce very small observable effects because they 
are hidden by the existence of unstabilities due to neutronization. 

8 Conclusions 

If we believe that the quantum vacuum in curved spacetime gives rise to some 
stress-energy, it is plausible that the gravitational effect of this contribution 
is equivalent to adding the term aR^ + hR^yR^" to the standard term, i?, 
of the Einstein-Hilbert action, at least if the curvature scalar R is not large. 
The presence of this term might give rise to effects observable, in principle, 
in terrestrial and the solar system observations. Then the present knowl- 
edge puts bounds on the possible values of the parameters a and 6, namely 
i/a, \/—h < 50 km. Also the vacuum stress-energy will produce some effects 
in the structure of Newtonian or slightly relativistic stars. In particular in 
white dwarfs it might produce stable hydrostatic equilibrium in stars above 
the Chandrasekhar limit, but the effect could not be observed because such 
stars would be unstable against neutronization. The effect of the vacuum 
stress-energy on the internal structure of the sun would be too small to be 
detected by observation. Our calculations suggest that the vacuum effects 
should be important in compact bodies like neutron stars. 
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